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We theoretically investigate the production of polarization-entangled photons through the biexci- 
ton cascade decay in a single semiconductor quantum dot. To accomplish a high degree of entangle- 
ment, despite the exciton finestructure splitting, one must either energetically align the two exciton 
states by means of external fields or erase the which-path information by post-selecting photons 
within the correct frequency range. Here we show that in the latter scheme, as well as in related 
proposals, the photon state becomes deteriorated through dephasing processes in the solid, and the 
degree of entanglement remains low despite filtering. 
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I. INTRODUCTION 

Conventional optics uses lenses and niters in order to 
manipulate light. The same elements are u sed for the 
manipulation of single or entangled photons.^ Indeed, 
in most cases of interest the working principle of passive 
optical elements is identical at the classical and quan- 
tum level, but there are situations were things are more 
cumbersome. Spectral filtering at the single photon level, 
for instance, affects the photon wavepacket and therefore 
modifies the quantum state of the photon. This has im- 
portant consequences for the creation of entangled pho- 
ton pairs in single semiconduct or q uantum dotsj 3 ^ which 
has recently been demonstrated^ and plays a key role 
in quantum information applications!^ 

In semiconductor quantum dots electrons and holes 
can be trapped to form excitonic complexes, which con- 
secutively decay through photon emission. Because of 
Coulomb correlation effects, the excitonic states usually 
have different energies and thus can be discriminated in 
optical experimentsP Quantum-dot based photon devices 
generally use th e last or last two photons emitted in such 
a cascade decayP^E^The last photon originates from the 
decay of the exciton, consisting of one electron and hole in 
the respective lowest unoccupied state, and the last but 
one photon from the decay of the biexciton state, consist- 
ing of two excitons with opposite spin orientations. The 
ambiguity of the biexciton decay, through either of the 
two ideally spin-degenerate exciton states, translates to 
a polarization entanglement of the emitted photons, that 
has been measured recently!^ 

The degree of entanglement is largely diminished by 
the electron-hole exchange interaction which splits the 
energies of the intermediate exciton states by a small 
amount and attaches a which-path information to the 
photon frequencies. Several experimental implementa- 
tions and proposals exist for overcoming this deficiency. 
First, the states can be br ought b ack to degeneracy by 
means of external magnetic 1 5 * 12 * 13 ! or electric^ fields, as 
well as through cavity modes with sufficiently strong cou- 
pling to both exciton and biexciton transitions.^ Alter- 
natively, spectral filtering can be used to post-select only 
those photons whose energy contains no which-path in- 



formation. 4 Finally, motivated by experimental observa- 
tions of quantum dots with a vanishing biexciton binding 
energy, it has been proposed to entangle photons pro- 
duced in different generations of the decay process .HSliZI 
As we will show in this paper, only the first protocol pro- 
vides a viable means for reaching a high degree of entan- 
glement. The other two inherently suffer from dephas- 
ing processes in the solid, which are always significant, 
and the degree of entanglement will consequently re- 
main rather low, thus limiting the performance of single- 
photon devices ESLL9| Qur results suggest that protocols 
for solid-state based quantum cryptography are more 
strict than previously thought. 

In this paper we provide analytical results for the de- 
gree of entanglement for the photons generated in the 
biexciton cascade decay in semiconductor quantum dots. 
We compare the different protocols based on spectral 
alignment of the intermediate exciton state, postselection 
through spectral filtering, and time reordering of photons 
produced in different generations. Our paper has been 
organized as follows. In Sec. |TT] we introduce our theo- 
retical model, and show how to compute the two-photon 
density matrix by means of the quantum regression the- 
orem. Section |III| presents the results obtained for the 
diffrcnt protocols and for realistic quantum dot parame- 



ters. Finally, in Sec. IV we draw some conclusions. 



II. THEORY 

A. Quantum state tomography 

In our theoretical approach, we consider the quantum 
dot level scheme depicted in Fig. la, consisting of the 
biexciton state u, the two exciton states with polariza- 
tions along a: and y, and the ground state g. By op- 
tical pumping j 4 * 5 * 10 * or electrical injection of carriers^^ 
the system is initially prepared in the biexciton state it. 
Thereafter it decays radiatively in a cascade process by 
emitting two photons. Because of the biexciton binding 
A, typically of the order of a few meV, the two photons 
have different energies and can be spectrally discrimi- 
nated. We shall refer to the photons from the biexci- 
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FIG. 1: (Color online) (a) Level scheme consisting of biexciton 
state it, exciton states polarized along x and y, and ground- 
state g. The biexciton energy is reduced by the biexciton 
binding A and the exciton states are energetically separated 
by the fmestructure splitting 8. H and V denote the polar- 
izations of the emitted photons, (b) Schematic of creation 
of entangled photons. The biexciton first decays by emit- 
ting a photon with polarization H or V, leaving behind a 
single exciton x or y in the dot, which decays in a second 
step. The ambiguity of the biexciton decay, through either 
of the two ideally spin-degenerate exciton states, translates 
to an entanglement of the emitted photonsP Because of the 
fmestructure splitting, the degree of entanglement is largely 
diminished. Several experimental implementations and pro- 
posals exist for overcoming this deficiency: (c) the states can 
be brought back to degeneracy by means of external fields^ 
(d) spectral filtering can be used to post-select only those 
photons whose energy contains no which-path information^ 
with h(u>) being the filter function, or (e) photons of different 
generations of the cascade can be entangled by bring ing th e 
biexciton transition in resonance with the exciton one! 16 * 17 ! 



ton and exciton decay as red and blue, respectively. If 
the intermediate exciton states are degenerate, the two 
decay paths are indistinguishable and the photons be- 
come entangled in polarizationP Such ideal performance 
is spoiled by the electron-hole exchange splitting^] as in- 
dicated in Fig.[l] To quantify the degree of entanglement, 
one usually performs a quantum state tomography^ by 
measuring the correlation between the two photons with 
a selected polarization combination, and averaging over 
the photon arrival times. This p roced ure provides us 
with the two-photon density matri x 21 ! 22 ! 



= avg 



I(t r )I(t b ) 



(1) 



which depends on the photon intensities I(t) of the first 
photon at time t r and the second photon at time t b . The 
elements of p^> 



(2) 



avg 



ElXt r )E- b (t b )E+(t b )E+Xt r ) 



then allow a complete reconstruction of the two-photon 
state. A and p, denote the horizontal and vertical photon 



polarizations H and V, respectively. Ef is the electric 
field operator for polarization A and for a positive or neg- 
ative frequency component.^ A schematic picture of this 
procedure is depicted in Fig. [2ja) for the diagonal ele- 
ments of the two-photon density matrix, and in Fig. |2[ 
for the off-diagonal elements. 

For the biexciton cascade decay s how n in Fig. 1 and in 
absence of exciton spin scatterings'^ the only non-zero 

(2) 

elements of the two-photon density matrix are p HH HH 

(2) 

and p HH yy, and those obtained by interchanging H with 
V. The entanglement of formation or concurrence,^ 
which provides a quantitative measure of the entangle- 
ment of the two emitted photons, is then given by the 
ratio of the off-diagonal and diagonal elements of p^ as 



C 



(2) 

Phh.w 



(2) 

Phh.hh 



(3) 



The concurrence is one for perfect entanglement and zero 
in case of no entanglement. 



1. Unfiltered case 

To establish a relation between the electric field oper- 
ators E^ 1 and the quantum dot dynamics, we employ the 
Green function G of classical Maxwell theory, which gives 
the field distribution for a localized current source. Let 
us consider the unfiltered case first. The expression^! 



Ef(t + t dclay ) = k z Gdf(t) 



(4) 

then connects the electric field with the dipole operator 
dr(i) of the excitonic transitions in the quantum dot. k 
is the wavenumber in vacuum. We have assumed that 
the photons arc emitted with a well-defined polarization 
A along the growth direction, as is th e case for quantum 
dots placed inside a microcavity.^E^ The main effect of 
the Green function is to introduce a time delay tdciay 
between the photon emission and its detection. We can 
now relate the photon-photon correlation function in Eq. 
([2]) to a dipole-dipole correlation function 



dx r {t r )d^(t b )d+(t b )d+ r (t r ) 



(5) 



where the dipole operators dj~ describe transitions be- 
tween different exciton states, associated with creation 
or destruction of a photon with given polarization A and 
color c. For simplicity, we have not indicated the con- 
stant time delay tdeiay between photon emission and de- 
tection. To arrive at the two-photon density matrix of 
Eq. ([2]), we have to multiply Q with the Green functions 
of Maxwell's theory, which describe the propagation of 
the photon fields from the dot to the photon detector, 
and to integrate over the photon detection times t r and 
t b . Note that the multiplication of Q with the Green func- 
tions G results in a mere constant, which drops out when 
calculating the two-photon density matrix or concurrence 
owing to the normalization condition tr[// 2 )] = 1. 
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2. Filtered case 

In the filtered case, the electric field is obtained from 
the expression 

£+(t + t delay )= / h{t-T)k 2 Gd+{T)dT, (6) 

Jo 

with h being the filter function. Upon Fourier trans- 
formation, the convolution becomes the familiar product 



h(u))d~^(u) of filter function and signal source. The terms 
under the integral describe how the photon wavepacket, 
created through excitonic decay, propagates away from 
the dot and becomes modulated by the filter [see Fig. 
|2jc)]. Due to the filter process, where the frequency in- 
formation is extracted from the temporal field evolution, 
p( 2 ) is determined by a correlation function, similar to 
Eq. jBJ, but now with all four dipole operators acting at 
different times 



h*(t r — T r )dT r / h* (tb — Tb) drb / h(tb — r' b ) dr' b I h(t r — t' t ) dr^. 



x (dx r (rr)d^(r b )d+ b (rl)d+ r (4) 



(7) 



For a broadband filter, whose spectral width is large in 
comparison to the excitonic lincwidth, the filter function 
becomes delta-like in time, and one recovers the dipole 
correlation function ^ of the unfiltered case. On the 
other hand, for a narrowband filter the integrals extend 
over a considerable time interval of the emission process. 



B. Quantum dot dynamics 

Next we show how to compute the dipole correlation 
function of Eqs. ^ and ([7]). The quantum dot dynamics 
is described as an open quantum system interacting with 
the environment. The free propagation of the quantum 
dot states is governed by the Hamiltonian 



H : 



E 

i=x,y 



Ei\i)(i\ + E u \u)(u\ 



(8) 



with i denoting the two different exciton states x and y 
with energy E x and E y7 respectively, and u the biexci- 
ton state with energy E u . Interaction with the environ- 
ment results in transitions between the different quantum 
dot states, and dephasing. We consider radiative decay 
and pure dephasing, which are expected to be the main 
scattering channels for excitons in quantum dots. Pure 
dephasing is due to phonon couplings^ and spectral dif- 
fusion due to charging centers in the vicinity of the dot. 
Other scattering channels, such as cross dephasin gj 13 ! 26 ! 
or exciton spin scatterings, are usually of minor impor- 
tance, and will be discussed at the end. We use the most 
simple description of radiative decay and pure dephasing 
in terms of scattering and dephasing rates "f r and jdi re- 
spectively. This will allow us to calculate the elements 
of p( 2 ^ analytically. Our rate equation approach is valid 
under the same conditions as the Born-Markov approxi- 
mation namely that the memory kernel describing the 
buildup of scatterings is sufficiently short in comparison 



to the scattering times. The memory of phonon dephas- 
ing, which is expected to b e the slowest process, is of the 
order of a few picosecond d 28 * 29 ! and thus much shorter 
than the dephasing time 7J 1 of several tens to hundreds 
of picoseconds. Quite generally, it is known that phonon 
dephasing cannot be fully described by a simple dephas- 
ing rate, and the lineshape of excitonic transitions con- 
sists of a no n-Lor entzian phonon background and a zero- 
phonon line! 25 * 28 ! The dephasing rate 7^ of our work mim- 
ics an averaged decay characteristics. We employ a mas- 
ter equation approach of Lindblad forrrPSl 



IP : 



where p is the density operator of the quantum dot states, 
and the Lindblad operators describe the various scat- 
tering channels listed in Table [T] For the Lindblad oper- 
ator associated with radiative decay, the last term in the 



TABLE I: List of Lindblad operators used in our calculations. 
7r and 7<j denote the radiative and dephasing rates, respec- 
tively. We assume that the biexciton rates are two times larger 
than the exciton ones, which is a reasonable approximation 
for both radiative decay and pure dephasingP^ The operators 
for cross dephasin g 13 * 2 * and exciton spin scatterings are given 
in the last two lines for completeness, but are not considered 
in the calculations given in Sec. [II] [except for Eq. £2}]. i 
labels the two exciton states x and y. 



Description 

Radiative exciton decay 
Radiative biexciton decay 
Exciton dephasing 
Biexciton dephasing 

Cross dephasing 
Exciton spin scattering 



Lindblad operator 

v^|0)(*l 
v^v |i)(w| 
JTd (|as><aj| + |i/><y|) 

V / ^7d| w )( w l 

\^nl*)(*l _ 
72 \y)(x\, v/72 \ x )(v\ 



Considered 

yes 
yes 
yes 
yes 

no 
no 
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FIG. 2: (Color online) Schematic of relation between quan- 
tum dot cascade decay and two-photon density matrix, (a) 

(2) 

Diagonal element of the two-photon density matrix p HH HH 
resulting from the decay of the biexciton u into the x polar- 
ized exciton and subsequently into the ground state g. The 
dipole operators d mediate the electric field amplitudes E, as 
described by Eq. We have indicated the states involved 
in the dipole transition rather than the polarization and color 
associated with the emitted photon. Wiggled lines indicate 
emitted photons, with polarization H or V depending on the 
orientation. Two electric field amplitudes build up a photon 
intensity Ihh at the detector (half circles). The two-photon 
density matrix is the correlation between measurements at the 
two detectors, as described in Eq. (fTj. The diagonal element 
Phh hh i s the result of averaging over the arrival times of 
the intensity Ihh of the first (red) and second (blue) photon, 
(b) Off-diagonal element of the two-photon density matrix 

(2) 

Phh w resulting from the decay of the biexciton u into a su- 
perposition of the x and y excitons. (c) p^ H w with a filter 
modifying the electric field amplitudes arriving at the detec- 
tors, according to Eq. |6|. (d) Time reordering by delaying 
the first generation of photons by a constant time, indicated 
by the thick black line between mirrors. 



parentheses of Eq. (|9| corresponds to the excitonic decay 
process, within which a photon is emitted from the dot. 

Introducing the shorthand notation Cp for the right- 
hand side of Eq. ([9]), where C is the Liouville superop- 
erator,^ we can write the formal solution of the master 
equation in the form 



p(t) = er* Ct Po 



{(^ £t )noo m +(^ £t ) Om }P0, (10) 



with po being the initial density operator. On the right- 
hand side of Eq. (10 1 we have split the time evolution 



ton emission, which changes the number of excitons. It 
will turn out that only the first part is needed in the cal- 
culation of the dipole correlation function ^ , which will 
allow us to considerably simplify our analysis. 

We next briefly discuss how to compute multi-time cor- 
relation functions of the form (Ai(ti)A 2 (t 2 ) • • ■ ), such as 
the dipole-dipolc correlation function, where A is an ar- 
bitrary operator. This is conve nient ly accomplished by 
the quantum regression theorem^^^ which allows to ex- 
press multi-time correlation functions in terms of single- 
time expectations. It implies that the fluctuations regress 
in time like the macroscopic averages, and holds exactly 
for systems initially decoupled from the environment and 
subject to a Markovian dynamics. Both approximations 
are valid for the description of the photon cascade decay 
in a quantum dot. For the initial decoupling it suffices 
to assume that the system will pass, at some point of the 
cascade, through the biexciton state, which consecutively 
decays through photon emissions. The applicability of 
the Markov approximation has been justified before in 
Sec. |IIB| According to the quantum regression theorem, 
a multi-time correlation function can be unraveled to a 
sequence of time evolutions, governed by e~ lCt , which are 
interrupted by abrupt transitions at times ti < t 2 < . . . 
according tcES 



(A 1 (t 1 )A 2 (t 2 ) . 



operator into two contributions, the first one associated 
with no photon emission, which does not change the num- 
ber of excitons, and the second one associated with pho- 



= tr ... A 2 e- lC{t ->- tl) Aie~ lCtl p 

Here Ap is a short-hand notation for Ap, if A is an oper- 
ator evolving with positive frequency such as d + , and for 
pA, if A is an operator evolving with negative frequency 
such as d~ . 

There is a subtle point regarding the use of the quan- 
tum regression theorem for the calculation of the dipole 
correlation function of Eqs. §5§ and ^ in case of a biexci- 
ton cascade, where the system is initially in the biexciton 
state po — \u)(u\. According to Eq. ( fil) , the dipole oper- 
ators d act from either the left or right hand side on the 
time-propagated pa, each of them creating a photon field. 
As we can create at most two photon fields from either 
side of po, it immediately follows that we have to con- 
sider in the time evolutions in Eq. (Ill only that part of 
e~' lCt associated with no photon emission. The remain- 
der corresponds to terms where more than four photon 
fields are involved, which are automatically zero in case 
of a pure cascade decay. 

All Lindblad operators considered in our calculations 
(see Table [l} either correspond to pure dephasing, which 
does not mix the different excitonic states, or to radiative 
decay. Thus, if the system is initially in the state |p)(<z|, 
with p and q labeling the quantum dot states 3, x, y, and 
u. the conditional time evolution (e~ lCt ) does not 

7 V / no em 

couple to any other state. From the master equation pM) 
we then obtain for the conditional time evolution within 
the usual rotating wave approximation^ 



5 



-iCt\ 



P Pq = [e 

f 1 f 

cxpl ift - 2tett 
expf ift - 2tett 
\ expf — iAi — 7tot* J 



|P)<9| = 

exp f — i^t — 
< p f — j 



cxp 
cxp I 



idt 
-i(A - 





(7r +7l)*) 







exp (-if t - 2tet t J cxp ( 

exp f — i5£ — (7 r -(- 71 ) f j cxp 

expf— 7 r ij cxp 

expf -i(A + 3 "2 Qt tl exp 



iAt - 7tot*] 
i(A- : 



i(A + 

-2 7 rt) 



I t- 



(12) 



where 7tot = 7r + Id is the sum of radiative and dephas- 
ing rates, 5 is the finestructure splitting, and A is the 
biexciton binding energy. The components of the condi- 
tional density matrix have one part associated with the 
free time evolution due to H, and another part associated 
with scatterings. The imaginary parts of the arguments 
of the exponentials are due to the different energies of 
the excitonic states, and the real parts due to scattering 
processes. Pure dephasing only affects those elements p pq 
which have an unequal numbe r of ex citons in p and q. We 
have included cross dephasing j 13 * 26 ! 71 in Eq. (12), which 
affects the elements p xy and p yx . This dephasing mech- 
anism will not be explicitly considered in our following 
calculations, but will be addressed in our final discussion. 



III. RESULTS 

A. Unfiltered case 

We now have all ingredients at hand to compute the 
elements of the two-photon density matrix. Let us con- 
sider the unfiltered case first. Following the prescription 
of the quantum regression theorem, Eq. (11), the corre- 



lation function is unraveled to a sequence of conditional 
density matrices p uu — ► p\^ — * p gg , which evolve under 
the condition that the system undergoes abrupt transi- 
tions, determined by the sequence of dipole operators, 
at times t r and tb- From the elements of the condi- 



tional density matrix of Eq. ( 12 ), we get for the diagonal 
and off-diagonal elements of the dipole-dipole correlation 
function 



Qhh,hh 
Qhh.vv 



We 
Ne 



-2 lr t r e -7r(*i,-tr) 

-2j r t r p -1r(,t b -t r )-i8(t b -t r ) 



(13) 

where Af is a normalization constant that includes dipole 
moments of different excitonic transitions. The entangle- 
ment of formation or concurrence [Eq. |3])] is obtained 
by averaging in Eq. ( 13 ) over the photon arrival times t r 
and tb, 



C 



It 



7,. + id 



(14) 



It has its maximum for zero finestructure splitting, and 
decreases with increasing 6, as shown in Fig. 3. This is 



because of the different phases acquired in the evolution 
of the intermediate exciton state, and the time average 
over the resulting phase factors. 
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FIG. 3: (Color online) (a) Concurrence as a function of 
finestructure splitting and for different dephasing rates, as 
computed from the expressions given in table [II] The solid 
line corresponds to the unfiltered case. The dashed (dotted) 
lines report results for 8ui = 1 peV (10 peV) which are plot- 
ted for positive (negative) finestructure splittings only. We 
use a radiative decay rate of f r — 1.6 peV which corresponds 
to a lifetime of about one nanosecond, (b) Concurrence as a 
function of filter width Su> and for S = 10 peV. In presence of 
dephasing the concurrence remains low despite filtering. The 
insets report density maps of the concurrence as a function of 
filter width and dephasing rate. 
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TABLE II: Elements of the two-photon density matrix and details of their calculation by means of the quantum regression 
theoremP^ In absence of exciton spin scatterings, the only nonzero elements are p^ H HH and p^ H yy, and those obtained by 
interchanging H and V. The concurrence^ is given by the modulus of the ratio between the off-diagonal and diagonal element!^ 
The correlation function is unraveled to a sequence of conditional density matrices, which evolve under the condition that the 
system undergoes abrupt transitions. Note that p' 2 ' is not p roperly normalized in case of filtering due to the loss of photons. 
In the entanglement of photons in different generations j 16 * 17 ! we assume that the photons originating from the biexciton decay 
are delayed by some time to. We use Ytot = 7r + 7d and c.c. for the complex conjugate of the preceding term. 

Figure Basis Ap Time order Decay path 

_ 

2 



Fig. la HH 0<t r <t b puu(t r ) — » p X x(t b — t r ) — » p gg 

HV < t r < t b Pira(tr) — * p X y(h — t r ) — * p gg 



7r 



2 7 r + iS 



Fig. Id HH a < T r < T r < n < r' b p U u(j r ) — > Pxu{t' t — r r ) — > p X x(j b - r r ) — ► p g x(r b — r b ) — > p gg 

< T r < T r < T b < T b p U u(T r ) ~* pux(T r - T r ) — * p X x(j b - T r ) — ► Pgx(r b — T b ) — > p gg 

< T r < Th < T r < t{, p U u(jr) — * pxu(j b — T r ) — * P g v.(j' r — T b ) — Y p g x(j' h — T r ) — > Pgg 

^ 1 ( 1 +^ + T , * ... | 4- ( .( 



8 Suj + iS + 37tot \5uj — i5 + "/tot Slj + iS + ytot {Sui + -ytot)(8uj — iS + y to t) 

HV" < T r < T r <T b <T b p U u(T r ) ~ * Pxv.{t' t - T r ) —> Pxy(j b — T r ) — > Pgy(T b — T b ) — > p gg 

< T T < T r < T b < T b puuij'r) —* puyijr — T r ) — > pxy(j b — T r ) — > P g y(j' b — T b ) — » p gg 
< T r < T b < t' t < T b puuijr) —* p X u(j b — Tr) —* Pgu(j'r — Tb) ~ > Pgy(j' b — Tr) — > pgg 

(Sij) 2 J r (21 



+ 



4 (Su + iS + 3jtot)(Suj + iS + Ttot) \7r + *<5 <5oj + 7 tot 
Fig. le HH 0<t r — t <t b p U u(t r — t ) p xx (t b — t r + t ) ^ p gg 

HV to < tr < tfc Puu(tr — to) —* pxu{t b — t r ) — > Pxy(t r ~ t b + t a ) — > P g y(tb — t r ) p gg 



to < tb < t r Puu(tb — to) — > puy(tr — tb) — > pxy(tb — tr + t() ) — > Px g (t r — tb) — > p 9£) 

7r ( o -1rt0 _ o -27tot«0^ 



7r + 27, 



"Other decay paths are obtained through r <-> r' and p pq — » p qp . 
For the off-diagonal elements one has to additionally replace a: <-» y. 



B. Filtered case 

The calculation of the two-photon density matrix in 
case of filtering is very similar but more cumbersome. 
We again employ the quantum regression theorem for its 
calculation, but we now obtain conditional density ma- 
trix elements p pq with unequal number of excitons, whose 
time evolutions arc affected by pure dephasing. To sim- 
plify our analysis, we assume that the filter function h(oj) 
has a Lorentzian shaped, with a central frequency posi- 
tioned in the middle of the two finestructure-split emis- 
sion peaks, and width Suj. Within this approach we can 
perform all integrals analytically. In the evaluation of 
the dipole correlation function |7]) we have to consider 
all possible time orderings for r, subject to the condi- 
tions t t < Tb and t' t < T' b which are a direct consequence 
of the cascade nature of the biexciton decay. Table |n] 
provides details about the resulting six time orderings, 
and lists the resulting diagonal and off-diagonal element 
of the two-photon density matrix (see also Fig. 2). 

Two limiting cases can be readily obtained from the 
elements of the filtered two-photon correlation function. 
First, for 6u> — > 00 the concurrence reduces to that of the 
unfrltered case. Second, for zero finestructure splitting 
the concurrence becomes one. This perfect entanglement 



is due to the absence of any which-path information in 
the photon frequencies. In the general case, the con- 
currence depends on the quantum dot parameters 5, 7 r , 
and 7d, as well as on the filter width Suj. Results for 
typical quantum dot and filter parameters are shown in 
Fig. 3. Panel (a) reports that in the filtered case the de- 
gree of entanglement first drops with increasing i5, and 
then increases again for larger values of S. This is due 
to the competition between the opposite trends of in- 
crease of which-path information with increasing 5, and 
the masking of which-path information through the fil- 
ter. The latter process is of importance when the two 
lines become energetically separated, but it comes at the 
price of a strongly reduced efficiency since in this regime 
practically all photons are absorbed by the filter. From 
the figure we observe that the concurrence dramatically 
drops in presence of dephasing. This is also evident from 
panel (b) which reports the influence of the filter width 
on the concurrence. For finite values of 7^ the degree of 
entanglement remains low despite filtering. This is due 
to the transition of the pure entangled photon state to 
a mixed, classically correlated photon state, which takes 
place in the process of filtering. 

The precise value of the dephasing rate depends on 
the detailed quantum dot parameters and on the cou- 
pling of the excitons to the solid state environment. The 
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line wid th of approximately 50 fieV measured in experi- 
ment!^ provides a hypothetical upper limit. This value is 
certainly too high, as it includes spectral diffusion which 
occurs on longer timcscales due to fluctuating environ- 
ment fields. These fluctuations affect the energies of the 
excitonic states but are not expected to have a drastic im- 
pact on the entanglement properties. For a lower bound 
of 7d we estimate a dephasing time of 340 ps, correspond- 
ing to a rate of approximately 5 /xeV, which wa s obtain ed 
in single-photon interferometer measurements.^ 19 * 33 ^ It 
is apparent from Fig. 3 that for such values of jd the 
concurrence is already restricted to small values. 



C. Time reordering 

Dephasing also plays an important role for entan- 
glement of photons created in different generations of 
the biexciton cascade (see Fig. le), which has recently 
been proposed for quant um dots with a vanishing biexci- 
ton binding energyASUU Using the notation of Fig. 1(e), 
within this approach one entangles the red photon fields 
H r and V r with the blue photon fields H b and V b . This is 
accomplished by delaying the fields H r and V b originat- 
ing from the biexciton decay by a constant time to [see 
Fig. [2^d)]. The elements of the dipole-dipole correlation 
function then read 



Ghhmh = (d Hr {t r - t )d Hh {t b )d^ h {t b )d^ r {t r - t Q ) / 
Qhh,w = (dH r {tr - to)dH b {tb)dy {t b - t )dy r (t r )ftl5) 



Table [TT| provides the results for the elements of the two- 
photon density matrix, which are obtained by evaluating 
Eq. ( 15 1 by means of the quantum regression theorem and 



averaging over the photon arrival times. Note that in this 
scheme the finestructure splitting d plays a different role 
than in the previous schemes, as it is used to distinguish 
the red from the blue photons. Employing again the ro- 
tating wave approximation, we consider this by setting 5 
equal to zero in the conditional density matrices ( 12 1. In 



experiments, S should be sufficiently large in order to al- 
low for an efficient filtering of the different photon colors. 



However, even in the case where filtering of the red and 
blue photon has no effect on the photon wavepacket, as 
assumed in the results given in the table, the time jitter of 
the photon emissions and dephasing significantly reduce 
the overall concurrence. Remarkably, even in absence 
of any dephasing the two photon wavepackets emitted in 
the cascade decay do not overlap completely, for any time 
delay to, and the concurrence is bound to values below 
one half. 



IV. CONCLUSIONS 

The radiative biexciton cascade decay in single semi- 
conductor quantum dots has the potential of being a 
source of triggered polarization-entangled photon pairs. 
Whether devices based on this mechanism will prove suc- 
cessful for quantum information applications depends on 
a number of issues. A high collection efficiency of pho- 
tons through the microcavity is crucial. In single-photon 
devices radiative scattering enhance ments by a factor of 
ten or more have been demonstratedJ^EH but things are 
more difficult in entangled-photon sources due to the dif- 
ferent photon colors. From a more fundamental perspec- 
tive, the issue of exciton finestructure splitting, which is 
a consequence of the general quantum principle of level 
repulsion, has to be addressed properly As we have 
shown here, through filtering^ or ent angle ment of pho- 
tons produced in different generation d 16 * 17 * the degree of 
entanglement becomes affected through dephasing losses. 
Such losses limit the device performance already at low 
temperatures, but make entanglement generation prob- 
ably impossible at elevated temperatures. In contrast, 
in the unfiltered case dephasing is a much less critical 
issue, with cross dephasing constituting the major cause 
of entanglement loss j 13 * 26 ^ and the devices could operate 
at higher temperatures.^ We thus believe that the tour 
de force approach of reducing the finestructure splitting, 
by means of growth optimization or external fields, is the 
most promising one for reaching a high degree of photon 
entanglement. 
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